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We investigate the initial progenitor mass dependence of the early-phase neutrino signal 
from supernovae taking neutrino oscillations into account. The early-phase analysis has 
advantages in that it is not affected by the time evolution of the density structure of the 
star due to shock propagation or whether the remnant is a neutron star or a black hole. 
The initial mass affects the evolution of the massive star and its presupernova struc- 
ture, which is important for two reasons when considering the neutrino signal. First, 
the density profile of the mantle affects the dynamics of neutrino oscillation in super- 
nova. Second, the final iron core structure determines the features of the neutrino burst, 
i.e., the luminosity and the average energy. We find that both effects are rather small. 
This is desirable when we try to extract information on neutrino parameters from future 
supernova-neutrino observations. Although the uncertainty due to the progenitor mass is 
not small for intermediate #13 (10~ 5 < sin 2 28n ^ 10~ 3 ), we can, nevertheless, determine 
the character of the mass hierarchy and whether #13 is very large or very small. 
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I. INTRODUCTION 

Core-collapse supernovae produce a huge flux of 
neutrinos of all types. The observation of neutrinos 
from SN1987A in the Large Magellanic Cloud was 
a milestone for neutrino astronomy [1,2]. Although 
some important constraints on neutrino properties 
were obtained [3-8], the total of 19 events in the 
Kamiokande and 1MB detectors was frustratingly 
sparse. On the other hand, SuperKamiokande is 
expected to detect about five to ten thousand neu- 
trinos from the next galactic supernova. Together 
with the other neutrino detectors arrayed around 
the world (e.g. SNO, LVD, ICARUS, IceCube), a 
high-statistics (nearby) supernova will provide an 
enormous amount of information on both super- 
nova and neutrino properties. 

Despite the new oscillation data from solar 
[9,10], atmospheric [11], reactor [12,13], and ac- 
celerator [14] neutrinos, the neutrino mixing an- 
gle #13 and the nature of the mass hierarchy re- 



main unknown. The constraints and implications 
of these data have been explored by many authors 
[15-23]. The greatest uncertainty for the study 
of neutrino oscillations in the supernova context 
is the spectral and temporal evolution of the neu- 
trino burst itself. Although the hierarchy of the 
average energies of the v e ,v e , and v x neutrinos 
{{E Ve ) < (Efj e ) < (E^)) is believed to hold, quan- 
titatively their specific spectra remain a matter of 
detailed calculation [24,25]. (y x denotes v^, v T and 
their antineutrinos.) Since a supernova is the end- 
product of a massive star, the properties of the 
neutrino burst will depend on the properties of the 
progenitor star, in particular its initial mass and 
envelope structure [26,27]. 

In this paper, we study how the initial mass of 
the progenitor star affects the early neutrino burst 
and the signature of neutrino oscillations in su- 
pernovae. The rest of the paper is organized as 
follows. In §11, we briefly review the evolution 
of massive stars, supernova explosions, and neu- 
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trino bursts. Conversion probabilities for progen- 
itor stars with various masses are calculated and 
compared in §111. Then, neutrino spectra with and 
without neutrino oscillations are shown in §IV. Fi- 
nally in §V we discuss the ramificatons of our re- 
sults and summarize our conclusions. 



II. EVOLUTION OF MASSIVE STARS AND 
SUPERNOVAE 

Here, we briefly summarize the evolution of mas- 
sive stars and supernovae. See [28] for a more ex- 
tended review. 

The evolution and death of a single massive star 
are determined by its initial mass and metallic- 
ity. For the metallicity of the Sun, a star with 
initial mass below about 8 Mq does not ignite 
carbon burning and forms a white dwarf. Above 
about 8 Mq, a star completes the advanced burn- 
ing stages, including for most stars silicon burning. 
This leads to the presupernova state, which is char- 
acterized by a core (iron or NeOMg) of roughly the 
Chandrasekhar mass surrounded by active burning 
shells and the accumulated ashes of oxygen, neon, 
carbon, and/or helium burning. The degenerate 
core in excess of the Chandrasekhar mass will dy- 
namically collapse and explode as a supernova. 

Up to an initial mass of 20 — 25M , the super- 
nova leaves behind a neutron star [29,30]. Almost 
all of the binding energy of the neutron star, 
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is radiated away as neutrinos. Here G, Mns and 
i?NS are the gravitational constant, the protoneu- 
tron star's mass, and its radius, respectively. Due 
to the differences in interaction strength, average 
neutrino energies vary with neutrino flavor. A 
full-scale numerical simulation by the Livermore 
group [31] using a multi-group flux-limiter diffu- 
sion approach for neutrino transport finds (E Ue } ~ 
13 MeV, (E Pe ) ~ 16 MeV, (E„ x ) ~ 23 MeV 
and almost perfect equipartition of the luminosi- 
ties [32] . Simulations with more sophisticated neu- 
trino transport [33], which are, however, limited to 
the early phase of neutrino burst, derive a smaller 
flavor-dependence for the average energy. These 
differences are essential if neutrino oscillations are 
to have a perceptible effect on supernova neutrino 
detection. (For details on supernova neutrinos, 



see, for example, the reviews by Burrows [34] and 
Suzuki [35].) 

For progenitor masses larger than 20 — 25M , 
sufficient mass may fall back onto the protoneu- 
tron star after explosion to turn it into a black 
hole. Even in this case, during the early phase 
which we consider in this paper, the properties of 
the neutrino burst are almost the same as stated 
above. 

The mass of the progenitor star affects the neu- 
trino oscillation signature of supernova neutrinos 
through differences in the mantle and core struc- 
tures. The density profile of the progenitor star, 
especially of the mantle, is important because it is 
related to the dynamics of neutrino flavor conver- 
sion. On the other hand, the structure of the iron 
core at the collapse determines the characteristics 
of the neutrino burst, e.g., the average energy and 
luminosity for each flavor. We discuss these in the 
following sections. 



III. NEUTRINO FLAVOR CONVERSION 
AND THE DENSITY PROFILE 

The evolution of a massive star is significantly 
affected by mass loss due to a stellar wind. In- 
deed, the mass loss can become so strong for a 
star with initial mass more than about 35 Mq and 
solar metallicity that the entire hydrogen envelope 
can be lost prior to the explosion of the star. It 
is suggested in [28] that the maximum in the fi- 
nal mass is about 20 Mq. In addition, the density 
profile of a star just before the supernova explo- 
sion will not depend much on the initial mass of 
the star. The final density profiles of stars with 
various initial masses are shown in Fig. 1 [36]. As 
is expected, they are similar just before collapse. 
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FIG. 1. Density profiles of stars just before super- 
nova explosion with initial mass 15Mq , 20M© , 40M© 
and 75M Q . 

The density structure in the range lOg/cc < p < 
10 4 g/cc is particularly important when consider- 
ing neutrino oscillations. Neutrinos produced in 
the high-density region of the iron core interact 
with matter before emerging from the supernova. 
Due to the non-zero masses and non-zero vacuum 
mixing angles among various neutrino flavors, fla- 
vor conversions can occur in supernovae. When 
the mixing angle is small, these conversions occur 
mainly in the resonance layer, where the density is 
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Am 2 is the mass squared difference, 9 is the mix- 
ing angle, E v is the neutrino energy, and Y e is the 
mean number of electrons per baryon. Since the in- 
ner supernova core is too dense to allow resonance 
conversion, we focus on two resonance points in 
the outer supernova envelope. One that occurs at 
higher density is called the H-resonance and the 
other, which occurs at lower density, is called the 
L-resonance. If the mass hierarchy is normal, both 
resonances occur in the neutrino sector. On the 
other hand, if the mass hierarchy is inverted, the 
H-resonance occurs in the antineutrino sector and 
the L-resonance occurs in the neutrino sector. 

The dynamics of conversion including large mix- 
ing case is determined by the adiabaticity param- 
eter 7, 
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which depends on the mixing angle and the mass- 
squared difference between the involved flavors. In 
eq. 3, n e is the electron number density. We define 
these as: 



#13 and Am 2 3 , at the H — resonance, (4) 



and 

#12 and Am 2 2 , at the L — resonance. (5) 

Here, the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix is taken as: 
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1,2,3 (i < j). When 7 » 1, the resonance is re- 
ferred to as an 'adiabatic resonance' and the fluxes 
of the two involved mass eigenstates are completely 
exchanged. However, when 7 <C 1, the resonance is 
called 'nonadiabatic' and the conversion does not 
occur. The dynamics of the resonance in super- 
novae is surveyed by Dighe and Smirnov [15]. 

The density profile of the star comes into 
the adiabaticity parameter 7 as the scale height 
n e /\dn e /dr\. A smaller scale height, that is, a 
steeper density profile results in less adiabatic res- 
onance. The scale heights of stars with various 
initial masses, calculated from the density profiles 
shown in Fig. 1, are given in Fig. 2 as functions 
of the density. The scale heights vary significantly, 
but differences between various initial masses are 
factors of 2 or 3 at the densities relevant for reso- 
nances. 
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FIG. 2. Scale heights n e /\dn e /dr\ of stars just 
before the supernova explosion with initial masses 
15 M ,2O M ,4O M Q , and 75 M Q . 

To estimate the effects of these differences on 
neutrino oscillations, we solve numerically the evo- 
lution equations for the neutrino wave functions 
along the density profiles shown in Fig. 1. We 
take the following values for the neutrino parame- 
ters: 

sin 2 20i2 = 0.84, Am 2 2 = 7 x l(T 5 eV 2 , 

= 3.2 x l(T 3 eV 2 . (6) 



sin 2 26» 23 = 1.0 
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The values of 6*12 and Amf 2 are taken from the 
global analysis of the solar neutrino observations 
and the KamLAND experiment [37] and corre- 
spond to the large mixing angle (LMA) solution 
of the solar neutrino problem, while those of #23 
,and Am 2 3 are taken from the analysis of the at- 
mospheric neutrino observations [11]. As for (9 13 , 
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we take sin 2 20 X3 = 1CT 4 . With this 13 , the H- 
resonance is neither perfectly adiabatic nor nona- 
diabatic and its adiabaticity is sensitive to the den- 
sity profile. The normal hierarchy is assumed, but 
the results are the same for the inverted hierarchy 
if v e is replaced by u e . 

In Fig. 3, the evolution of the probability 
P(y e — > v e ) that a neutrino emitted as a v e at the 
neutrinosphere remains a v e is shown. The neu- 
trino energies on the plot are 5 MeV and 40 MeV. 
The H-resonance radii and the final probabilities 
can be quite different for different progenitors and 
neutrino energies. The observationally important 
quantity is the final conversion probability and this 
is shown in Fig. 4 as a function of the neutrino 
energy. The differences are not so small, about 
0.05 — 0.1 at all energies. But if sin 2 26\% is very 
large or very small, the adiabaticity parameter is 
very large or very small, respectively. In this case 
such difference of scale height as in Fig. 2 will 
not affect the neutrino conversion probability. The 
range of sin 2 26* 13 is shown to be sin 2 26* 13 < 10~ 5 
and sin 2 26>i 3 > 10" 3 in [23]. 
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FIG. 3. Evolution of probabilities P(v e — > v e ) that 
a neutrino emitted as a v e at the neutrinosphere re- 
mains a «/ e at some radius. Upper and lower curves 
correspond to neutrino energies of 40 MeV and 5 MeV, 
respectively. 
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FIG. 4. Energy dependence of 

probabilities P(^ e — > v E ) that a neutrino emitted as 
a v e at the neutrinosphere remains a v e at the surface 
of the star. Also shown are the perfectly adiabatic and 
non-adiabatic cases for the H-resonance. 

Although the density structure evolves as the 
shock wave propagates, it takes about 2 seconds 
for the shock wave to reach the H-resonance re- 
gion [39] . We concentrate our analysis on the early 
phase when the shock propagation effect can be 
neglected. The potential time dependence of neu- 
trino oscillations due to shock propagation is dis- 
cussed by several authors [38,39,17,40]. 



IV. NEUTRINO BURSTS AND NEUTRINO 
OSCILLATIONS 

As was stated in the previous section, massive 
stars experience significant mass loss. For cur- 
rent empirical mass loss rates, all solar-metallicity 
stars initially more massive than about 35 Mq 
arc thought to become hydrogen-free objects of 
roughly 5-M© at the end of their thermonuclear 
evolution. The corresponding upper limit to the 
mass of the final iron core is about 2M© [28]. 

The mass of the iron core is determined roughly 
by the Chandrasekhar mass, which for a zero- 
temperature, constant Y e , Newtonian structure is 
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However, there are numerous corrections, some of 
which are large [41]. To a first approximation, the 
non-zero entropy of the core is important and 
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is the electronic entropy per baryon. More massive 
stars have higher entropy and contain larger iron 
cores on average. However, this general tendency 
is moderated by the loss and redistribution of en- 
tropy that occurs during the late burning stages. 
Thus, the mass of the iron core as a function of the 
initial mass will be somewhat uncertain in that a 
small change in the initial mass results in a large 
difference in the iron core mass. According to [28], 
the mass of the iron core is 1.2(1.4) — 1.6M when 
the initial mass is beween 10(20)M Q and 40M Q . 
This weak dependence of the iron core mass on 
the ZAMS progenitor mass leads to a somewhat 
universal neutrino burst. 

Figs. 5-10 show the evolution of the average neu- 
trino energy and number luminosity in the early 
phase up to 200 milliseconds after bounce. The 
calculation is based on dynamical models of core- 
collapse supernovae in one spatial dimension, em- 
ploying a Boltzmann neutrino radiation transport 
algorithm, coupled to Newtonian Lagrangean hy- 
drodynamics and a consistent high-density nuclear 
equation of state. Details of these simulations are 
described in [33]. As can be seen, the major fea- 
tures of the early neutrino burst are almost inde- 
pendent of the initial mass. 
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FIG. 6. Evolution of the number flux at the Earth 
of v e neutrinos due to a neutrino burst from a progen- 
itor of initial mass 15M© , 20M Q , 40M© and 75M© at 
a distance of 10 kiloparsecs (kpc) . Time at the bounce 
is set to zero. 
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FIG. 5. Evolution of the average v e energy of 
a neutrino burst from a progenitor of initial mass 
15M© , 2OM , 4OM and 75M Q . Time at the bounce 
is set to zero. 
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FIG. 7. Same as Fig. 5, but for v £ 
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FIG. 8. Same as Fig. 6, but for i/ e . 
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FIG. 9. Same as Fig. 5, but for u x 
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FIG. 10. Same as Fig. 6, but for v x . 

Combined with the discussion in the previous 
section, we conclude that the mantle structure 
and the features of the neutrino burst depend lit- 
tle on the initial mass of the progenitor star if 
sin 2 26»i3 % 10~ 5 or sin 2 20 13 > 1(T 3 . On one 
hand, this means that we can not easily obtain 
information about the initial mass from observa- 
tions of neutrinos during the first 200 milliseconds 
after bounce. On the other hand, this situation 
is desirable for extracting information about the 
neutrino parameters themselves. 
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FIG. 11. Time-integrated event spectra at Su- 
perKamiokande using our supernova model with initial 
mass 2QMq at a distance of 10 kpc. 
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FIG. 12. Time evolution of the event number in Su- 
perKamiokande at a distance of 10 kpc using our su- 
pernova model with an initial mass of 20Mq . 
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FIG. 13. Time-integrated event spectra at SNO and 
a distance of 10 kpc using our supernova model with an 
initial mass of 20Mq . Only the charged-current events 
due to the process u e + d^p + p + e~ have been taken 
into account. 
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FIG. 14. Time evolution of the event number at 
SNO and 10 kpc using our supernova model with an 
initial mass of 2QMq . Only the charged-current events 
due to v e (y e + d — > p + p + e~) absorption have been 
taken into account. 

In our previous papers [20-23], we investigated 
the possibility of extracting information on neu- 
trino parameters, especially f3 13 and the neutrino 
mass hierarchy, from future supernova neutrino ob- 
servations. We used a numerical supernova model 
from the Livermore group [31,32] with an initial 
mass of 15Af . Since the average neutrino energy 
is flavor dependent (E Ve < Ey c < E^), neutrino 
oscillations make the spectra of observed v e and v e 
neutrinos harder. The extent to which observable 
neutrinos become harder depends on the adiabatic- 
ity parameters at both the H- and L-resonances, 
which depend on neutrino parameters. We con- 
sidered models with four sets of neutrino param- 
eters: normal-LMA-L, normal-LMA-S, inverted- 
LMA-L and inverted-LMA-S. Here, normal-LMA- 
L means that the mass hierarchy is normal (mf < 
mf. < m§) and 9 13 is large (sin 2 29 13 = 0.043), 
inverted-LMA-S means that the mass hierarchy g=" 
is inverted (to 2 <C to 2 < m|) and #13 is small 
(sin 2 20i3 = 10~ 6 ), and so on. The other parame- 
ters are the same as in the previous section. 

In Figs. 11-14, we show for a 20- M Q progenitor 
the time-integrated event spectra and the time evo- 
lution of the event number at SuperKamiokande 
(SK) and SNO. The events at SNO include only the 
charged-current events due to the process v e + d — ► 
p+p + e~ , while all relevant processes are included 
in the SK event estimates. The detector properties 
and cross sections used to calculate these events 
are described in our previous paper [21]. It should 
be noted that the neutronization burst events at 
SNO are significantly suppressed due to neutrino 
oscillation. 

In [21,23], we used as a measure of neutrino 
oscillation the ratio of high-energy event number 
to low-energy event number at SuperKamiokande, 
whose events are dominated by v e , and SNO, which 
can identify v e . 

High (low) means higher (lower) than 20 MeV. 
Our conclusions were that inverted-LMA-L is 
clearly distinguishable from the other models al- 
though the distinction between normal-LMA-L 
and normal- and inverted-LMA-S is rather diffi- 
cult. In addition, we determined that normal- 
LMA-S and inverted-LMA-S are completely degen- 
erate [23] . One of the uncertainties in our analysis 
was that we did not previously take into account 
the initial progenitor mass dependence. 



We plot in Fig. 15 the above-mentioned event- 
number ratios for SK (-Rsk) and SNO (i?SNo)- In- 
cluded are results for the above four neutrino pa- 
rameter models, as well as for the no oscillation 
case, and various initial progenitor masses. Since 
we use only early-phase events (< 0.2 sec), the 
statistical errors are large, whereas in our previous 
papers we incorporated the full evolution of the 
neutrino burst based on the simulations of the Liv- 
ermore group [31,32,26,27]. It can be seen that dif- 
ferences due to initial mass are much smaller than 
those due to neutrino parameter and statistical er- 
rors. Thus, we can conclude that although the un- 
certainty due to the unknown progenitor mass is 
not small for intermediate #13 (10~ 5 < sin 2 2# 13 < 
10~ 3 ), we can nevertheless draw conclusions about 
whether #13 is very large or very small and about 
the mass hierarchy. 
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FIG. 15. Plots of Rsk and Rsno for the four neu- 
trino parameter models mentioned in the text, as well 
as for the no oscillation case. Results for various ini- 
tial progenitor masses are given. The error bars show 
statistical errors only. 



V. DISCUSSION AND SUMMARY 

In this paper, we have investigated the initial 
mass dependence of the early-phase neutrino sig- 
nal from supernovae. The early-phase analysis has 
the advantage that it is not affected by the time 
evolution of the density structure of the star due to 
shock propagation and is independent of whether 
the remnant is eventually a neutron star or a black 
hole. The initial mass does not affect neutrino os- 
cillations in two senses: the density profile of the 
mantle, which is important for neutrino flavor con- 
version, and the final iron core structure, which 
determines the features of the neutrino burst, are 
almost independent of the initial mass. This is de- 
sirable when we try to extract information from fu- 
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turc supernova-neutrino observations on neutrino 
parameters. Although the uncertainty due to the 
progenitor mass is not small for intermediate #13 
(1CT 5 < sin 2 2#i3 < 10~ 3 ), we can nevertheless ex- 
tract information about whether #13 is very large 
or very small and about the neutrino mass hierar- 
chy. 
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